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ABSTRACT 


Modifications were made to Corgnati's (1971) steady state, symmet- 
rical hurricane model, which was based on a model proposed by Riehl 
(1963), and in which bulk aerodynamic transfer formulas were used to 
predict the air-sea interaction below a specified vortex. The modifi- 
cations resulted in a more flexible model in which hurricane intensity 
was controlled by sensible and latent heat fluxes, which were in turn 
dependent on the ocean surface temperature. A boundary layer model 
based on similarity theory was then used, to specify air-sea interaction, 
by replacing the bulk transfer formulas with Cardone's (1969) extension 
of Blackadar's (1965) two-layer neutral, baroclinic boundary layer 
model. To achieve realistic results the time scale was modified to 
include centripetal acceleration, which resulted in a decreasing height 
of frictional influence with decreasing radius. Varying Cardone's for- 
mulation of roughness length and the ratio of eddy exchange coefficients 
for heat and moisture resulted in stress and heat flux values in agree- 
ment with values determined from hurricane observations. Including the 
complete boundary layer model resulted in greater sensitivity to ocean 
temperature variations than when bulk transfer formulas were used to 


determine the air-sea interaction. 
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I. INTRODUCTION 


It is well recognized that sea surface temperature affects hurricane 
intensity; the hurricane winds in turn, exert a stress on the ocean sur- 
face and induce the extraction of heat and moisture. A "Steady state" 
hurricane model requires that a balance exists in this energy exchange 
between air and ocean. 

Riehl (1963) described a two-layer model for an axisymmetric hurri- 
cane in which a balance existed between the heat source and momentum 
sink at the air-sea interface. The "steady state" upper level vortex 
was then specified based on the air-sea interaction within the boundary 
layer. 

Corgnati (1971) used Riehl's model to construct a numerical model 
which used predictive equations for heat and moisture within the 
boundary layer, and bulk transfer formulas for fluxes at the air-sea 
interface. It was felt that considerable improvement in Corgnati's 
model could be attained through the use of recent boundary layer models 
based on Monin-Obuhkov "similarity theory"; specifically through the 
use of Cardone's (1969) extension of Blackadar's (1965) two-layer neutral 
baroclinic boundary layer model. 

The purpose of this study was to modify Corgnati's numerical model 
so as to better define the air-sea interaction. 

Several numerical simulations of hurricane development have been 
made with complicated multi-level models. These models include some 
interaction with the ocean. The scope of this study was to model only 
the boundary layer beneath an upper level vortex, defined by Riehl's 


model, which responds to the air-sea interaction. 
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II. DISCUSSION OF THE MODEL 


A. REVIEW OF RIEHL'S MODEL 

Riehl's two-layer model consisted of an inflow layer and an outflow 
layer (Figure 1). Riehl and Malkus (1961) have shown the conservation 
of absolute angular momentum is a reasonable first approximation in the 
outflow layer. This constraint leads to the expression for computing 
the radius at which the tangential component of the outflow goes from 
cyclonic to anticyclonic, 


= (1) 


= Jl 
R, = (2 a R, /£) 


where Ri is the eye wall radius and uS is the tangential wind at Ri. 
i 
Riehl and Malkus (1961) also demonstrated the approximate conserva- 


tion of potential vorticity in the inflow layer. Riehl (1963) applied 


this to the equation of motion in polar coordinates (see Figure 2) 


dV VV OT ’ 
e er - 3 Ber 1 e 
TOME A p R90 " p oz" 32) 


and integrated over the depth of the inflow layer to get 


Bis = constant , (3) 


where = is the tangential surface stress. Using the expression 


= 2 
T, = CP. I bé (4) 


and assuming the drag coefficient (Cp) and surface density (6,) to be 
constant, also that the inflow angle is small enough so that cosine A 
can be closely approximated by unit yields the expression 


u Ri /2 = constant . (5) 
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Figure 1. Illustrating hurricane mass circulation V_ is tangential 
wind component, 9. is equivalent potential temperature, R is radius. 





Figure 2. Illustrating polar coordinate system used, Va is tangential 
wind component, V, is radial wind component, R is radius, 8 is angular 
measurement. 
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Riehl (1963) and Shea (1972) have shown this to be a reasonably good 
approximation to the VTR profiles. 


Equation (2) can also be solved for the radial wind profile, 


2 
V c O 1 (6) 
r AP V JA 2 

u 


where G is the acceleration of gravity, and f is the coriolis parameter. 

Integration of the hydrostatic equation, for typical values of 
equivalent potential temperature c in a hurricane, produces the 
relationship 


P' = - 2.56 on > (7) 


where P' is the pressure departure from 1005 mb and o is the departure 
from 350K. Riehl (1963) also determined a relation between the maximum 
tangential wind Ci and the equivalent potential temperature (8) at 
the eye wall by ee Ye to the height gradient of the E 
surfaces and applying a um relationship to obtain 


Vo, (m/sec) = 14.1(8, - 8 ee (8) 
1 1 O 


where subscript "o" denoted the value at Ro: 

From these equations it can be seen that the upper level vortex is 
dependent on the increase in En in the boundary layer. The problem 
becomes one of modeling the boundary layer air-sea interaction so that 


the A increase corresponds to the upper level vortex. 
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B. CORGNATI'S MODEL 
1. Prediction Equation 
a. Atmospheric Prediction 

In the atmosphere the predicted variables were heat and 
moisture in the form of potential temperature (8 ) and specific humid- 
ity (q,) respectively. These two parameters, along with pressure (P), 
determined 9 which was the one factor in determining the maximum wind 
(Equation 8). 

Partial derivatives were approximated using finite differ- 


encing techniques to solve the thermodynamic equation in the following 


form SS RO. ae „Ss G 
Tes O C AP , (9) 


where w is the vertical motion and Qs is sensible heat transfer. The 
first term on the right hand side was neglected due to symmetry. The 
boundary layer was assumed well mixed with height, and therefore, 9 
could be calculated from the surface values of air temperature and 
pressure. Moisture changes were determined in a Enti fashion using 


the expression oq, OV ola oV Rd, Ea Q G 
2|+ — , (10) 
"SP DA 


=. ne ROR 


where Qo is the latent heat transfer. 
Corgnati (1971) used bulk aerodynamic transfer formulas for 
determining the transfer of latent and sensible heat in the form 
Q, 7 $6, 6,6, - TOV , (11) 
Qe » 3m ne 5 TA 4 2) 
where VD is the specific heat of dry air, L, is the latent heat of 
vaporization of water, Py is the atmospheric density at the sea surface, 


I, - T, is the air-sea temperature difference, qu and q, are the 
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specific humidity of the air and the saturated specific humidity at the 
sea surface temperature. The surface wind (V) was approximated at 70% 
of the wind at the top of the boundary layer (Equation 5). Furthermore, 
Robinson's (1966) heat transfer coefficient (2/3 Cy) was used, which 
takes molecular processes near the surface into account. 

Vertical motion at the top of the boundary layer was com- 
puted from the continuity equation in polar coordinates, 


OV ORV ST 


e r 


R58 t BR + T =0 (13) 
where the first term is zero due to symmetry. Vertical motion at the 
surface was assumed to be zero, and a boundary layer 100 mb thick was 
used, based on the assumption that effects within the boundary layer 
would not exceed this height. 

b. Oceanic Prediction 

Predicted variables in the oceanic portion of Corgnati's 
model were sea surface temperature (2) and the saturated specific 
humidity at the sea surface temperature. A temperature gradient of 2C 
per 100 meters with a layer depth of 100 meters described the ocean 
structure. The amount of heat removed, through latent and sensible 
heat transfer, was computed and the temperature of the ocean surface 
reduced accordingly.  Corgnati found that, since the model was 
stationary, reductions in T, led to reduced air-sea transfers and, 
therefore, reduced wind speed. In order to test only the atmospheric 
portion of the model the sea temperature was fixed. 


2. Initialization of Surface Temperatures and Humidities 


The relationship between 0. , tangential wind at the eye wall 
i 
(V, ), and pressure at the eye wall (P1) made it possible, by specifying 
T 
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NE , to solve for a and Pi using Equations (7) and (8). The surface 
i i 
pressure distribution was then determined by assuming a gradient wind 


balance. 

Initialization of the temperature and moisture fields, both 
in the boundary layer and at the ocean surface, were ee parts 
of the model. The surface water temperature was set at 30C, then the 


available moisture at the ocean surface was calculated from *he inte- 


grated form of Clapeyron's equation 





e 
S s T-273 T 


where e, is the saturated vapor pressure. Air in contact with the 
water was assumed Saturated, with the specific humidity calculated from 


4 


-622 e 
S 


do" P- .373e, ° (15) 


The equivalent potential temperature was computed by applying the 
relationship between height gradient of the isobaric surfaces and the 
tangential wind field. 

By specifying the initial state from the equivalent potential 
temperature at the eye wall, the effect of expansion cooling and 
sensible heat source at the surface were taken into account. Riehl 
(1963) suggested that e must rise above 350-352K to allow deep con- 
vection, which warms the column and lowers the central pressure to 
values which can sustain hurricane force winds. 

Specific humidity of the air (q,) was calculated assuming the 


initial air temperature (T,) equal to the water temperature, from 


C INS 
4, = T, T Ln T >? (16) 
v a 


where LS is the pseudo-equivalent air temperature. 
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C. MODIFICATIONS TO CORGNATI'S NUMERICAL MODEL 


Corgnati (1971) used a 17 by 17 rectangular grid with 15 nautical 
mile spacing. The center of the grid was the hurricane center, and 
variables were computed along a diagonal and interpolated to the other 
grid points. Due to the symmetry of the model, the four quadrants had 
identical results, and as much information could be derived through the 
use of a one-dimensional model. Resolution was increased by using 100 
grid points spaced at 3 kilometer intervals along a radial. One dis- 
advantage of the increased resolution was that smaller time steps had 
to be used, although this was partially compensated by the reduction in 
the number of grid points. Inherent instability in a forward-in-tima, 
centered-in-space finite differencing technique, made it necessary to 
change to a leapfrog technique with a forward time step every 24 time 
steps. 

The equivalent potential temperature gradient was the determining 
factor in predicting the maximum wind. When deriving Equation (8), 
Riehl assumed cyclostrophic flow, so that the relationship was only a 
good approximation in the area of maximum winds. Assuming gradient 


winds led to the relationship. 


2 
A ae E oes 
R e p OR ° 


< 


(17) 


which was used during the initialization to determine the equivalent 
potential temperature field. 

The interdependence of the equivalent potential temperature field 
and the wind field makes the selection of the eye wall and outer radii 
important factors in determining the wind field. The inner radius was 


allowed to adjust towards the radius of maximum equivalent potential 
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temperature, rather than remain fixed as in Corgnati's model, and the 
outer radius was determined from Equation (1). The maximum tangential 
wind component was then determined from Equation (17) using the gradient 
of equivalent potential temperature from the eye wall to the outer 
radius. 

Corgnati (1971) found that relative humidity values at the eye wall 
region exceeded saturated conditions after a short period of time. He 
theorized that cumulus clouds, in the form of "hot towers", with the 
moist ascending air being replaced by relatively dry subsiding air, 
would serve to remove this excessive moisture. To model the "hot towers" 
Corgnati removed all moisture in excess of a critical value of 95% rela- 
tive humidity, and measured this as convective precipitation. A similar 
problem was encountered at the outer radius, where the computed value of 
o was found to be significantly higher than the 359K specified by Riehl. 
Since the maximum wind was dependent on the radial gradient of AE the 
increased 9 caused a decrease in the wind field. To correct this, a 
critical relative humidity of 75$ was used at the outer radius. This 
simulated the relatively dry subsiding air outside a hurricane. Several 
profiles of critical relative humidity were experimented with and will 
be discussed later. 

Corgnati (1971) experimented with the release of latent heat, within 
the boundary layer, of the excess moisture removed through cumulus con- 
vection. Releasing all the latent heat within the boundary layer greatly 
increased the hurricane intensity. The bases of cumulus clouds within 
hurricanes have been observed to be lower towards the eye. To model 


this, latent heat release was made an inverse function of radial distance. 


18 





With these changes it was now possible, through latent heat re- 
lease, relative humidity control and ocean temperature adjustments, 


to model steady state hurricanes of different intensities. 


D. CARDONE'S BOUNDARY LAYER MODEL 

The use of bulk aerodynanic transfer formulas requires some know- 
ledge of heat and moisture exchange coefficients as well as winds 
within the surface layer. Recent advances in the study of atmospheric 
turbulence have led to a greater confidence in describing the distri- 
bution of wind and temperature in the atmospheric boundary layer. This 
progress has arisen largely through the application of the Monin-Obuhkov 
"similarity theory", which has been supported by observational studies. 
Cardone (1969) has extended Blackadar's (1965) two-layer, neutral, 
baroclinic boundary layer model, over fixed terrain, to the boundary 
layer over the sea. 'McConathy (1972) evaluated Cardone's formulation 
of the marine boundary layer. The results of McConathy's study were 
used to apply Cardone's specification of the boundary layer, instead of 
the bulk aerodynamic transfer formulas in Corgnati's model, and to 
determine a radial wind profile throughout the boundary layer. 

The Monin-Obuhkov "similarity theory" for the constant stress, or 
surface layer, predicts that a universal relation should exist for the 


wind shear, 


EN E (18) 


and the temperature gradient, 


98, 9r 
SEU Ne (19) 
where the friction en (U,) is defined as 
1/2 
2 (20) 
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and the scaling temperature 
1 H-: 


Cp = 
* Der 


(21) 


k is von Karman's constant of .35, z is the height, H is the total heat 
flux, b, and b. are functions of stability. Matching the integrated equa- 
tions for the Ekman layer to the integrated equation for the constant 
stress layer at top of the surface layer led to the following set of 
equations, in terms of U, and inflow angle (a), 


U 


E" [2kB sina 0, à? , (22) 
and 
U,  ky2 sin G - a) 
2 See a S (23) 


de LnB Ro - HEL) 


where the integrated profile stability parameter is defined by 








: u -1 T a ^ 
(21-6 - 3Lnp + 2Ln( A ES 
and the height of the surface layer (h) is defined by 
BV 
Hes cee, | (25) 
the surface Rossby number (Lettau 1959) by 
V 
Ro = A ; (26) 


and the modified stability length (Panofsky 1963) by 


K, JU Qu /dz 


Ds ME 


K kcoü/z  ^* (27) 


In these expressions = is the surface geostrophic wind, and By is 
Blackadar's constant of 0.0003. The drag coefficient (C) was defined 


in terms of the roughness parameter (Z,) and stability parameter as 
2 


= k 
vg T i Zu rcm = 
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The roughness parameter was formulated by Cardone (1969) as 


_ 6.84 x 107? 


o U, 


4 


Z + 4.28 x 10°? Ue - 4.43 x 10" (29) 


where the second term on the right hand side represents Charnock's(1955) 
formilation for aerodynamically rough flow under high wind speed condi- 
tions, 2 
.035U, 


Zo = eo > (30) 


Equations (22) and (23) were solved simultaneously for UN and A from 
the external parameters s f, and air-sea virtual temperature differ- 


ence. 


E. APPLICATION OF CARDONE'S BOUNDARY LAYER MODEL 

Cardone's model was applied to the balanced hurricane to determine 
the fluxes of heat and moisture as well as momentum transfer. The 
determination of a mean radial wind component, using Cardone's model, 
assured continuity in determining advection within the boundary layer. 

The tangential wind, which represented a gradient wind at the top 
of the boundary layer, was corrected for thermal wind and used as the 
surface gradient wind. 

To include the centripetal acceleration, the coriolis parameter was 
replaced by 


: (31) 


By assuming the turbulent transfer coefficient for heat (K, ) equals the 
turbulent transfer coefficient for momentum (Ka? and knowing L'from 


the solution of Equations (22) and (23), the total heat transfer was 


3 


determined using Ka UC p T 
L's =L s-——t 43 , (32) 
Ka kGH 
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The total heat (H) was then separated into latent and sensible components 


by assuming their ratio was determined by 


Q CNO e B 
es =) (33) 
da 7 8, 


The predictive Equations (9) and (10) required a radial wind for 


Qe Ly 
advecting the heat and moisture. The assumption that heat and moisture 
are homogeneous, within the depth of the boundary layer, made it possible 
to use an integrated value of radial wind. This was determined by 


solving the integrated form of Equation (18) 


cd 


2o 2 e 
U k [L nZ WE DU 9 (34) 


ZI 


at z = i=1, 2, 3, 4, and averaging these values. The mean radial 
wind in the surface layer was then determined by multiplying by the sine 
of the inflow angle (X). The inflow within the spiral layer was deter- 


mined through integration of the Ekman solution through the spiral layer 


to get co 
se 3T 3T : 
\ V „dz= /2 e sinX z——- [cos(C4 7 -Bh)-sin(a+"> -Bh)], (35) 
h dn 4 4 
where 
5 0172 
B= GU kh] > 


The two integrated values of un from Equations (34) and (35) were then 
used in determining the vertical motion (Equation 13), and the advection 
of heat and moisture (Equations 9 and 10) in the predictive portion of 


the model. 





III. DISCUSSION OF RESULTS 


The numerical process encompassed three stages. First was an ini- 
tialization phase in which the model was started with a given Vo? T 2 
T; and Ris then the model went through an adjustment process in which 
the moisture content in the boundary layer increased until either a 
balance was reached between latent heat flux and moisture advection, or 
the moisture reached the critical relative humidity values. The final 
stage was a "steady state" in which the time variations of heat and 
moisture in the boundary layer were small. 

The following discussion deals with the interpretation of results 
obtained when various parameters were changed, both in the initial and 
adjustment stages of the model. Due to the complex interactions of the 
model, a simplistic approach of establishing a base run and then varying 
one parameter at a time was used. This type of experimentation made the 


effects of the changes more obvious. 


A.  CORGNATI'S MODIFIED MODEL 

The modifications to Corgnati's model, as discussed in the previous 
section, were incorporated to achieve more flexibility in describing 
hurricanes of varying intensity.  Corgnati discussed in detail the adjust- 
ment process. He found that steady state was dependent on latent heat 
release and critical relative humidity used to simulate the cumulus con- 
vection. These parameters control, to a good part, the heat and 
moisture at the eye wall (Ri); which in turn determines the equivalent 
potential temperature (ou and, therefore, the hurricane intensity 


dl 
(Equation 8). 
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The parameters varied were initial air temperature, water temperature, 
critical relative humidity, and latent heat release within the boundary 
layer. Figures 3 through 15 illustrate the relationships between the 
various parameters, and were obtained from the base run with water tem- 
perature set at 30C, air temperature initialized at 29C, critical rela- 
tive humidity set at 95% at R; > with a linear profile to 75% at Ro? and 
no latent heat release within the boundary layer. 

Air temperature (Figure 3) exhibited a nearly isothermal (less than 
.3C) profile fron R, to about SR,» with an approximately logarithmic 
drop from SR, to a minimum temperature at Ri- The temperature drop cor- 
responds to the very strong horizontal pressure gradient (Figure 4). 

The adiabatic cooling due to expansion would cause about a 5 C drop, 

but, as Riehl (1954) pointed out, the blowing spray in the high wind 
region would greatly increase the effective surface area and, therefore, 
the latent and sensible heat transfer. This enhanced transfer would off- 
set the adiabatic cooling and account for the slight temperature decrease 
towards the eye as observed by Deppermann (1939). The bulk ee 
formulas (Equations 11 and 12) do not take into account the enhanced 
transfer due to the blowing spray. However, the increased wind speed 

and air-sea temperature difference were sufficient to increase the 
sensible and latent heat fluxes and offset much of the adiabatic cooling. 

Relative humidity (Figure 5) was at the critical value which indi- 
cates the steady state stage was one in which the flux and advection 
of moisture was kept in balance through the removal of moisture exceed- 
ing the critical relative humidity. Figure 15 indicates the moisture 


removed to keep the relative humidity at the critical value. 
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Figure 3. Radial profile of air temperature for the base run 
of Corgnati's modified model, T. - 30C, T, initialized at 29C. 
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Figure 4, Same as Figure 3 except pressure. 
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Figure 5. Same as Figure 3 except relative humidity. 





Figure 6. Same as Figure 3. except specific humidity. 
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Specific humidity (Figure 6) followed the relative humidity in the 
isothermal region. The sharp decrease in temperature reduced the 
amount of moisture, despite the increase in critical relative humidity, 
and the result was a peak in specific humidity outside the eye wall. 

Equivalent potential temperature (Figure 7) is dependent on tempera- 
ture, moisture, and pressure. In the nearly isothermal region and with 
small pressure gradient, A is governed by the amount of moisture. From 
3R, to Ri the temperature decrease has a two-fold effect in that it also 
reduces specific humidity and 9, becomes linear. The gradient of 8 
between Ri and Ro determines e (Equation 17), and Equation (5) de- 

i 
scribes the tangential wind profile (Figure 3). Pressure and wind are 
linked through the gradient wind equation, and the sharp wind increase 
is reflected by the strong pressure gradient. 

In a hurricane the winds rotating about the eye experience both 
centripetal and coriolis acceleration (Figure 9). In the outer portion 
of the hurricane coriolis dominates. However, as V increases and 
radius decreases, centrifugal force becomes dominant, supporting the 
strong pressure gradient in the eye wall regions, and causing a decrease 
in integrated s The decrease in M in the near surface region requires 
cross isobaric flow; this results in a frictional force opposing the 
pressure force. The importance of these various forces and their region 
of influence will become more apparent in the discussion of results 
obtained using Cardone's boundary layer model. Figure 10 shows the 
effect of cyclostrophic flow on the radial wind. The strong vertical 
motioa (Figure 11) at the eye wall is due to the decrease in radial wind 
towards the eye wall and the decreasing operates (Equation 13). 


In reality the vertical motion is due to deep convection maintained 
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Figure 7. Same as Figure 3 except equivalent potential temperature. 
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Figure 8. Same as Figure 3 except tangential component of the wind. 
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Figure 9. Illustrating the forces acting in the boundary layer 
of an axisymmetric hurricane. 





Figure 10. Same as Figure 3 except radial component of the wind. 
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Same as Figure 3 except vertical motion. 
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through the buoyant energy gain represented by the 2 increase from R, 
to Ri. The convection is sufficient to cause the large pressure gradi- 
ent which is maintained by the cyclostrophic balance. 

Latent heat flux (Figure 12) is dependent on wind speed and the 
difference in specific humidity of the air and saturated air at the sea- 
surface temperature (Equation 12). The lov moisture content of the air 
at R5 causes a local maximum in heat flux. The increase in moisture 
dominates the region from R, to SR, causing a decrease in flux despite 
the increasing wind. Maximum heat flux occurs at R; where the decrease 
in pressure and temperature acts to increase the specific humidity 
difference and the wind reaches a maximum velocity. 

Sensible heat flux (Figure 13), which is dependent on wind speed 
and air-sea temperature difference (Equation 11), reaches a maximum at 
Ri and is about an order of magnitude less than the latent heat flux. 

Precipitation was simulated by two mechanisms., One was the large 
scale moisture flux through the top of the boundary layer (Figure 14), 
which reflects the specific humidity and vertical motion. The second 
mechanism is the removal. of moisture which exceeds the critical relative 
humidity (Figure 15). Both curves are accumulated during the first 15 
hours and, therefore, are influenced by the adjustment phase. The peak 
in Figure 15 at Ro is due to the movement of RS during adjustment. 

l. Experiments with Temperature Variations 

In a series of experiments the horizontally uniform water tem- 
perature was varied from 28 to 31C while air temperature was 
initialized at 28C. The radial profiles of the various parameters were 
similar to those obtained from the base run and illustrated in Figures 3 


to 15. There were, however, changes iu magnitude and significant values 
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Figure 12. Same as Figure 3 except latent heat flux. 





Figure 13. Same as Figure 3 except sensible heat flux. 
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Figure 14. Same as Figure 3 except total moisture flux through the 
top of the boundary layer for first 15 hours. 
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Figure 15. Same as Figure 3 except for total moisture removed by 
critical relative humidity mechanism for first 15 hours. 
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are presented in tabular form (Table 1). Air temperature (T, at Ro 
o 
adjusted to within one degree of the water temperature, and the air 
temperature (T, ) at R; was about 1.5C cooler than T - One of the 
o 

interesting relationships was the dependence of the model on the differ- 
ence in e Eon RS to Ri. Even when the water temperature was fairly 
cool, the increase in zs due to air-sea exchange was large enough to 
sustain hurricane force winds. Since 9 is representative of the ambi- 
ent values in a tropical atmosphere, nds of 350K are necessary 
for deep convection, pe was limited to a minimum value of 350K for 


O 
determining va (Equation 8). This limitation is reflected in the value 


of Ye for the 28C water temperature case. One of the assumptions of the 


i 
model is that the energy source for the hurricane is the ocean; there- 


fore, the hurricane should have more energy available and increase in 
intensity over warmer water. Intensity did increase from 49 m ae for 
the 28C ocean to 62 m for the 31C ocean. 

2. Experiments with Variations of Critical Relative Humidity 

Critical relative humidity was used to simulate decreases of 

moisture in the boundary layer due to convective mechanisms. To test 
the effect of the critical relative humidity, three profiles were speci- 
fied (Figure 16). Resulting relative humidity profiles (Figure 17) show 
that in the convex case relative humidity did not reach critical values 
except in the eye wall and outer regions. In the concave and linear 
cases relative humidity reached critical values at all radii between Ro 
and R; Latent heat flux is dependent on the difference in vapor 
pressure of saturated air at sea surface temperature and ambient air; 
because of this dependence the latent heat flux distribution was affected 


by the specification of critical relative humidity (Figure 18). 





Table I 


Significant values in the modified model when the 
water temperature was varied from 28 to 31C, and 
air temperature was initialized at 28C. 


T T T 8 8 V P 

W a. a es E e. O 

Oo o- o ol o = 
coa) ¿Em (K) (K) (m/sec) (mbs) 
28 26.7 21:9 359 347 49 975 
29 27.0 28.5 363 350 57 965 
30 27.6 20.2 367 353 59 961 
31 28.4 29.9 371 356 62 958 
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Figure 16. Radial profiles of the three specified critical 
relative humidity curves. 
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Figure 17. Radial profiles of relative humidity resulting from 
the three different specifications of critical relative humidity 
Ty = 30C, T initialized at 29C. 
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Figure 18. Same as Figure 17 except latent heat flux. 
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Figure 19. Same as Figure 17 except equivalent potential temperature. 


37 


Hurricane intensity, which is determined from 2 at Ro and Rj» did not 
vary. The eye wall was made to adjust towards the maximum 95 in the 
convex case maximum ds occurred outside the eye wall and Ri adjusted to 
36 km. In the concave and linear cases maximum 9 occurred at or inside 
the eye wall and R; adjusted to 24 km (Figure 19). 

It was concluded that the actual form of the critical relative 
humidity curve between fixed extrema was not crucial in predictii the 
intensity of che storm; however, it does affect the distribution of 
heat exchange. 

3. Experiments with Variations of Latent Heat Release 

Assuming dry adiabatic expansion in the inflow would result in 
air temperature in the eye wall region as much as 5C colder than ob- 
served temperatures. Riehl (1954) hypothesized that latent heat release 
or sensible heat flux, enhanced by blowing spray, must offset the expan- 
sion cooling. To simulate the effects of this enhanced heating, the 
latent heat of condensation of the moisture released to maintain relative 
humidity below the critical value (Figure 15), was added to the boundary 
layer. Since the enhanced heating would be a function of the amount of 
blowing spray, and the amount of blowing spray would be dependent on 
wind speed, and wind speed is a function of radius, the percentage of 
heat released was made a linear function of radius with the maximum value 
at Ri; and no latent heat release at Roe Corgnati experimented with re- 
leasing 100% of this latent heat at all radii and found a marked increase 
in intensity. Table II illustrates the effects of latent heat release in 
the modified model. The increase in 9,» while 9 remained unaffected, 


i o 
resulted in an increased intensity. 
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Variations of significant values with latent 


Table II 


heat released in the boundary layer. 


latent heat 
release at Ri 


(%) 


T 
a 


T 
o 


Ge) 
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8 


e. 
1 


8 V P 
e e, O 
o 
(°K) (m/sec) (mb) 
354 58 963 
354 69 . 951 
354 70 944 
354 72 941 
354 74 938 


Although latent heat was used, the heat gain could also be pro- 
duced by an enhanced sensible heat transfer due to blowing spray. In 
all probability both latent heat release and enhanced sensible heat 
transfer play some role in offsetting expansion cooling in the boundary 
layer. The cr profile of the model, without latent heat re- 
lease, was in agreement with observed temperature profiles. The model 


was therefore run with no latent heat release in the boundary layer. 


B. CARDONE'S BOUNDARY LAYER MODEL 

Cardone's two layer baroclinic boundary layer model assumes hori- 
zontal homogeneity, stationarity, Ka equal to a constant in the Ekman 
layer and 


R (37) 


in the surface layer. The model also involves various empirically 
determined constants based on data from mid-latitudes. Despite these 
limitations it was felt that a boundary layer based on similarity theory 
would provide an internally consistent model of an hurricane boundary 
layer. The inputs to Cardone's model were surface gradient wind, 
virtual temperature difference between the air in the boundary layer and 
saturated air at the sea surface temperature, and latitude. 
1. Application of Cardone's Boundary Layer Model 

Several problems arose with the direct application of Cardone's 
model in place of the bulk transfer equations and the integrated inflow 
from the tangential wind equation. The model hurricane decayed rapidly 
during the adjustment phase, and heat and moisture fluxes were well 
below the 3000 to 5000 cal a m values determined by Whitaker 


(1967) to occur in the hurricane force wind region of hurricane Betsy. 
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Investigation showed that radial wind (Figure 20), determined from Equa- 
tions (34) and (35), increased inward from 3R, to Ri. Because M. must 
go to zero in the eye, vertical motion at the eye wall was on the order 
of 1500 mb n In this region the vertical transport of moisture ex- 
ceeded the sum of the horizontal transport and evaporation, which de- 
creased a and led to a decrease in intensity.  Haltiner and Martin 
(1957) EN shown that the balance of forces in a vortex may be depend- 
ent on centrifugal force. This was included in Cardone's model by 
replacing f by f' as formulated in Equation (31). The time scaling 


factor then became 


(f ee : (38) 


which had the effect of reducing the height of frictional influence as 
va increased. In particular, the height of the constant stress layer 
was reduced from about 50 meters at R, to about 20 meters at R; (Figure 
21). To maintain consistency within the boundary layer, the centripetal 
acceleration also had to be included in term B (Equation 36) of the Ekman 
spiral. The formulation for V, within the spiral layer became 

VS 172 AC 1/2 


f+ f+ 
Í R e : 3T R 
hU, k2 Z oe pis n U, k? Z | (39) 








V, (z)- Jle 


where V, (2) is defined by 





Pe CO 1/2 V, (2) 1/2 

E f+ —— 

V (z)= V TJ2 e R Z sin Q co Q+ 3T c - 

o e h U, k2 z 4\hu, k2 /? (40) 


To determine ue in the spiral layer an iterative technique was employed 
in determining V, (2) using heights (z) at 10 meter intervals. Then V (z) 
r 


was determined at each height (Figure 22) using the Vo Le) for that height. 
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Figure 20. Radial profile of radial component of the wind resulting 
from direct application of Cardone's boundary layer model. 
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Figure 21. Radial profile of height of the surface layer, illustrat- 


ing decreasing height of frictional influence due to cyclostrophic 
flow. 
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Figure 22. Vertical variation of radial component of the wind in 
the spiral layer, derived using iterative technique. 
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The integrated radial wind in the spiral layer was then determined by 
averaging the values of v,(z) over the spiral layer. Using an iterative 
technique proved to be rather time consuming. Because V, (2) increased 
rapidly with height, the assumption that V, 62) was equal to Vo at the top 
of the boundary layer resulted in integrated values of V. (z) that closely 
approximated the averaged values of v (z) (Table 3). The assumption 

that Vo was constant with height in the spiral layer was used in deter- 
mining the integrated value of vr (Figure 23). Including centrifugal 
forces was now consistent within the entire model and the rapid decay 
described earlier, as well as the extreme vertical motion, was no longer 
evident. 

In a study of water temperatures in the wake of hurricane Betsy, 
Whitaker (1967) estimated the heat flux in the hurricane force wind 
region to be between .035 and .058 cal ens e The application of 
Cardone's model produced heat fluxes with a maximum value around .02 cal 
En? E (Figure 24). One of the assumptions in Cardone's model was 
that the ratio of the turbulent transfer coefficient for heat (Kp? to 
turbulent transfer coefficient for momentum (K) was unity.  Businger et 
al (1971) found "The ratio of the eddy diffusivities (a, ) for heat and 
moisture is greater than unity. The value at neutral is 1.35; and shows 
a marked increase with instability." Since heat transfer, as determined 
in the model, is directly dependent on an (Equation 32), increasing a, 
increased the heat flux. In a hurricane the large quantity of blowing 
spray should serve to enhance the heat transfer. Furthermore, the 
velocity imparted to this blowing spray should serve to increase the 


momentum transfer. To simulate the increasing heat flux and momentum 


transfer, the roughness length (Z,) was increased by increasing the 


44 





Table III 


Comparison of integrated V_ values obtained 
using iterative technique, and assuming Vo (Z) 
constant. 


Radius (km) 30 60 90 120 150 180 210 


V_ (m/sec) 
iterative 3.14 3,2273,42 3,49 3,40 3.38 3,20 
technique 


V (m/sec) 


3,22 3.21 3.47 3753 3743 3.401372 
NO const 


180 
R km 





Figure 23. Radial profile of integrated radial component of the 
wind assuming cyclostrophic flow and ns constant within the spiral 
layer. 
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Figure 24. Radial profile of combined sensible and latent heat flux 
obtained using Cardone's model, C, = .00428, = 1,0, T_ = 30C, 

2 ah w 
and T, initialized at 29C. 


an 






4 
Q,+aQ, 
x 10? 


cal/c msec 


R. 60 120 180 240 R 300 
: R km 


Figure 25. Same as Figure 24 except a, = 240, C} = .05. 
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the coefficient (C5) of the us term in Cardone's formulation of Z , (Equa- 
tion 29). It was found that realistic values of heat flux were obtained 


when C, was increased from .00428 to .05, and a 


9 taken as 2.0 (Figure 25). 


h 


Cardone's model, with cyclostrophic flow, C, increased to .05, 


2 
and an increased to 2, was now used to determine the sensible and latent 
heat fluxes, replacing Equations (11) and (12), and the average radial 
wind component, replacing Equation (6). 

A base run was made in which water temperature was set at 30C, 
air temperature initialized at 29C, critical relative humidity set at 
95% at R; with a linear decrease to 75% at Ro» and no latent heat re- 
leased within the boundary layer. This base run was compared to the 
base run illustrated by Figures 3 through 15 and discussed earlier.  Car- 
done's model produced greater latent and sensible heat flux (Figures 12 
and 25) which led to slightly warmer temperatures at Ri and increased 
9. As a result of the increased 9L? od was about 5 m eu greater. 
Ele 4 compares some of the ee le of the base runs. The 
use of Cardone's boundary layer model provided some additional para- 
meters that were indicative of the structure of the boundary layer, as 
illustrated in Figures 26 through 29. The modified stability length 
(L'), illustrated in Figure 26, can be interpreted to be proportional to 
the ratio of mechanical and buoyant production of turbulent kinetic 
energy. The large increase of L', with increasing Nas reflects an in- 
crease in the mechanical versus the buoyant production of kinetic energy. 
Surface stress values (Figure 27), which increases rather dramatically to 
a maximum value at Ri, are in close agreement with stress values deduced 


by Miller (1962) using data from hurricane Helene (Table 5). The increase 


in stress towards the center is due to the increase in wind and the 
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Table IV 


Comparison of significant values obtained from 
two base runs. 


Corgnati's Cardone's 
parameter modified model model 
Y (m/sec) 58 62 
i 

i Ce) 27.9 28.0 
T 

T mo 29.5 29.2 
P o 

8 - 8 ( K) 13 15 
e. e 
i o 

En (mbs) 963 953 
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Figure 26. Same as Figure 25 except modified stability length. 
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Figure 27. Same as Figure 25 except surface stress. 
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Table V 


Comparison of surface stress determined 
by Miller, and values obtained from the 
hurricane model. 


surface stress Radius 0-20nm  20-40nm 40-60nm 
qae Miller's 60.0 5150 27.6 
values 
Radius 20nm 30nm 50nm 
computed 60 35 20 
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rougher surface. The increase of stress with decreasing radial distance 
is reflected in the increase in inflow angle within the constant stress 
layer (Figure 28); also in the radial velocity at the top of the sur- 
face layer (Figure 29). These values are comparable to inflow angles 
and radial velocities observed in actual hurricanes. One of the assump- 
tions made in Corgnati's model was that the cosine of the inflow angle 
could be closely approximated by unity. Despite the increase in the 
inflow angle, this appears to be a good approximation. 
2. Temperature Experiments with Cardone's Model 

To further test the use of Cardone's boundary layer model, ex- 
periments were conducted in which water temperature was varied from 28 
to 31C. These experiments were similar to the ones conducted with 
Corgnati's modified model, and discussed in section IIIAl. The heat 
flux distribution obtained, using Cardone's model, appeared to be more 
sensitive to temperature en than the heat flux distribution 
obtained using the bulk aerodynamic transfer formulas (Figure 30). Due 
to the differing heat flux, the variation in the difference between 
extrema values of ER was greater, using Cardone's model, than the 
variations obtained using bulk aerodynamic transfer formulas. The 
result of these greater variations was that, using Cardone's model, a 
sea surface temperature increase from 28 to 31C produced a In increase 
OF 25 m E CMS (Table 6); while using bulk aerodynamic — formulas, 
the same temperature increase produced only a 13 m s increase in 


Vo . It was felt that the greater sensitivity to sea surface tempera- 
i 


ture, exhibited by Cardone's boundary layer model, more closely simulated 


the dependence of a hurricane on the oceanic heat source. 
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Figure 28. Same as Figure 25 except inflow angle in surface layer. 
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Figure 29. Same as Figure 25 except radial component of the wind 
at top of surface layer. 
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Figure 30. Radial profiles of combined latent and sensible heat 
flux, cedo values obtained from Cardone's model with T = 31C 
and T, = 29C, and Corgnati's modified model with T, = 31C and 
I7 W290, T, initialized at 28C. i 


Table VI 
Changes in enews values due to vary- 


ing water temperature, air temperature 
initialized at 28C. 


Ww ay a ce e e. o 


Co Co CO er eec) OD) 


23 2/7.) 2580 358 348 45 979 
29 12522 02932 353 349 58 960 
30 27.8 28.7 368 351 66 952 
3105092875 12953 70373 359 70 944 
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IV. CONCLUSIONS 


The results of this study have indicated that a marine boundary 
layer based on similarity theory can be applied to an hurricane and 
improve the effects of air-sea interaction. 

First, results indicate that ane time scaling factor in cyclo- 
strophic flow is dependent on (f + A rather than f as in 
Cardone's model. The effect of the centripetal acceleration was to 
reduce the height of frictional influence as wind speed increased and 
radial distance decreased. 

Second, observed temperature fields indicate that enhanced heat 
flux is necessary to offset the adiabatic expansion cooling within 
the boundary layer of an hurricane. This enhanced heat flux is likely 
to be due to the large quantities of blowing spray observed in an 
hurricane. Consequently, the ratio of the coefficients of heat dif- 
fusivity to momentum diffusivity was increased from unity to 2. 
Increasing the roughness length, as a function of friction velocity 
squared, was used to simulate the increased momentum transfer and 
resulted in surface stress values that appear to be realistic. The 
increase in roughness length was further justified by resultant heat 
flux values that were in agreement with determined values. 

Third, the removal of moisture from the boundary layer, through 
convective mechanisms, appears to have a strong influence on air-sea 
interaction. Vertical advection by the mean motion is not sufficient 
to export the large horizontal flux of moisture within the boundary 


layer and the evaporation. The removal of moisture in excess of a 


54 


critical relative humidity value appears to adequately simulate "hot 
tower'" convection and surrounding subsidence. 

Finally, in the earlier model of Corgnati (1971) it kad been neces- 
sary to specify a drag coefficient ad extrapolate a surface wind to 
determine air-sea interaction. However, the external parameters 
required to use similarity theory were surface gradient wind, and the 
virtual temperature difference between saturated air at the sea surface 
temperature, and the air in the boundary layer. The vertical structure 
of the wind and an implied drag coefficient were inherent in the 
iterative solution based on similarity theory. 

The more complete boundary layer specification based on similarity 
theory appears to be justified in view of the good agreement obtained 
between predicted and observed or estimated values of radial inflow, 
heat flux, and surface stress. Furthermore, the boundary layer model 
calculates wind profiles within the surface layer, as well as surface 
stress; and effects of both stability and changes of surface roughness 
on air-sea interaction have been included in the hurricane model. 

One of the assumptions made in this study was that the stability 
length was defined by density stratification based on virtual tempera- 
ture. This assumption was crucial in the determination of latent heat 
flux. 

Further studies, in the area of air-sea interaction under high wind 
speed conditions, are necessary to provide guidance in the selection of 
empirical constants in the similarity solution. Specifically, the 
effects of blowing spray, in enhancing heat flux and momentum transfer, 


may be important in maintaining tropical cyclones. 
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